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We study the time-dependent dynamical properties of two-component ultracold fermions in 
a one-dimensional optical superlattice by applying the adaptive time-dependent density matrix 
renormalization group to a repulsive Hubbard model with an alternating superlattice potential. 
We clarify how the time evolution of local quantities occurs when the superlattice potential 
is suddenly changed to a normal one. For a Mott-type insulating state at quarter filling, the 
time evolution exhibits a profile similar to that expected for bosonic atoms, where correlation 
effects are less important. On the other hand, for a band- type insulating state at half filling, the 
strong repulsive interaction induces an unusual pairing of fermions, resulting in some striking 
properties in time evolution, such as a paired fermion co-tunneling process and the suppression 
of local spin moments. We further address the effect of a confining potential, which causes 
spatial confinement of the paired fermions. 
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Ultracold atoms in optical lattices offer unprecedented 
potential for exploring many-body effects. 1 This system 
is considered to be an ideal realization of the Hubbard 
model that has been intensively studied for many years in 
condensed matter physics. Recent successful demonstra- 
tions of the Mott-insulator phase transition 2-4 prove that 
we can obtain a comprehensive understanding of compli- 
cated many-body effects from a quantitative comparison 
of experiments and theories. 

High controllability of experimental parameters is one 
of the outstanding features of ultracold atoms. In op- 
tical lattices, this allows us to study the nonadiabatic 
dynamics of ultracold atoms by measuring the time evo- 
lution of atomic states after a quick modification of the 
lattice potentials or a sudden change in inter-atomic in- 
teractions via Fcshbach resonance. Experiments are car- 
ried out in an ultrahigh vacuum so that energy is con- 
served and the quantum dynamics can remain coher- 
ent during the time evolution. Various intriguing prop- 
erties have already been reported: the collapse and re- 
vival of a Bose-Einstein condensate, 5 the quantum ver- 
sion of Newton's cradle, 6 the co-tunneling of bound atom 
pairs, spin dynamics mediated by exchange interac- 
tions, 8 and second-order atom tunneling. 9 These exper- 
imental observations have stimulated theoretical inves- 
tigations. Several numerical approaches have been ap- 
plied to the Bose-Hubbrad model 10-14 and the Fermi- 
Hubbard model 15 ' 16 to reveal the nonequilibrium dy- 
namics of many-body states that is largely unavailable 
in traditional condensed matter systems. However, most 
theoretical investigations have focused on bosonic sys- 
tems owing to the current experimental situation. This 
naturally motivates us to study the nonadiabatic dynam- 
ics of fermions with strong correlations. 

In this paper, we investigate the nonadiabatic dy- 

*E-mail address: yamamoto@tp.ap.eng.osaka-u.ac.jp 



namical properties of two types of insulating states in 
a one-dimensional (ID) optical superlattice with two- 
component correlated fermions. First, we show the time 
evolution of the insulating states when there is a sud- 
den change in a superlattice potential. We clarify how 
the particle correlations emerge in the nonadiabatic pro- 
cess. In particular, we clarify that the unusual pairing 
of fermions caused by the strong repulsive interaction 
produces a characteristic time-dependent profile of the 
physical quantities. The effect of a harmonic confining 
potential is also addressed. 

Let us start with an extended Hubbard model, which 
may describe two-component ultracold fermions in a ID 
optical superlattice with two-site periodic potential, 



n = -J E ( c !«rC*+V + H - c ) + ^E 



-v d 



E ■ 

i—odd.a 



L + l 



(1) 



where c\ a {ci^ a ) is a creation (annihilation) operator of a 
fermion at the i-th site with spin a{=], J.). Here J is the 
nearest-neighbor hopping integral, L the total number of 
lattice sites, U (> 0) is the on-site repulsive interaction 
and Vd is the strength of the two-site periodic superlat- 
tice potential. V c is the curvature of a harmonic confin- 
ing potential, which will be dealt with later in this pa- 
per. The ground state properties of the system without 
a confining potential were clarified in connection with 
quasi-lD organic molecules, 17-20 and the effect of con- 
finement was addressed in the context of cold fermions. 21 
Among several intriguing phases caused by the superlat- 
tice structure, 17-21 we here consider two types of insulat- 
ing states as an initial state of time evolution: Mott-type 
and band-type insulating states, which are respectively 
realized at quarter- and half-filling with a large potential 
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(a) t=0 (V d /J>0) (b)t>0(V d /J=0) 
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Fig. 1. (Color online) Schematic diagrams of nonadiabatic con- 
trol of one-dimensional optical superlattice: (a) an initial two-site 
periodic superlattice with a large potential difference V d and (b) 
a normal lattice after a sudden disappearance of Vd for t >0. 



difference, V d ^> J- 

In the following, we deal with the dynamics of the 
above fermionic superlattice at quarter- and half-filling 
(we temporarily set V c = 0). To realize nonadiabatic con- 
trol of our system, we suddenly change the lattice poten- 
tial from V d /J = 20 (t = 0) to V d = 0(t > 0) as shown in 
Fig. 1. We make use of the adaptive time-dependent den- 
sity matrix renormalization group (adaptive t-DMRG) 22 
method which is an extended version of the original 
DMRG 23 supplemented by the time-evolving block deci- 
mation algorithm. 24 The adaptive t-DRMG method can 
deal with the real-time dynamics of strongly correlated 
systems. We thus calculate the time-dependence of the 
density profile (n,), the variance of local spins AS* 
(={S* )-(Sf) 2 ) and the double occupancy of fermions 
(TH,t n i,i) - We check the accuracy of t-DMRG by com- 
paring the results with the time-dependent exact diago- 
nalization for small systems. In the following, we adopt 
J as units of energy and time is measured in the unit of 

l/h. 

First, we investigate the nonadiabatic time evolution 
of a Mott insulating state at quarter-filling. Fermionic 
atoms are initially located at every other site because 
of the large potential difference Vd(=20) between even 
and odd sites, which may be described by a sequence of 
empty and singly occupied sites 0,1,0,1, . . .0. Thus the 
system is a variant of a Mott insulator realized in a su- 
perlattice. 17-21 In Fig. 2, we show the time evolution of 
local quantities after V d has been suddenly turned off. 
In Fig. 2(a) (Fig. 2(b)) we show the time-dependent local 
density (n,) at a given site i = 18 (t = 19) that is unoccu- 
pied (singly occupied) initially at t = 0. We find that the 
timc-cvolution profile is little affected by the interaction, 
and its oscillation period and amplitude are similar to 
those observed for bosonic atoms (approximately given 
by the Bessel function of the first kind). 11,12 Note that 
(rij) at even and odd sites always satisfies the conser- 
vation law of total fermions: the sum of (rij) in the unit 
cell gives unity. Although particle correlations do not ap- 
pear to be very important for (rij), they are somewhat 
relevant as regards the variance of local spins ASf and 
the double occupancy (rij^rij^). Figure 2 indeed shows 
that the double occupancy is suppressed and its oscil- 
lation period becomes shorter as U increases. Moreover, 
the double occupancy evolves differently at odd and even 
sites (Fig. 2(a) and (b)). Let us first look at the even 
site in Fig. 2(a). For a small U, fermions move rather 
freely, so that (rij^ri^i) is roughly given by the square 
of (rij). However, for a larger U, the system needs the 
additional energy U to realize double occupancy at even 




Fig. 2. (Color online) Time dependence of the local quantities for 
several choices of on-site interaction for the system at quarter 
filling with L = 36 and N = 18. Plots of local density (nj), 
variance of local spins AS 1 ^, and double occupancy (rij f rij i ) 
from top to bottom. The on-site interaction is chosen as U = 1 
(solid line), 5 (dashed line), and 9 (dotted line) at (a) i = 18 and 
(b) i = 19. 



sites, which thus reduces its probability and shortens the 
oscillation period (see discussions below for the oscilla- 
tion period). For odd sites, as seen in Fig. 2(b), the am- 
plitude of (ni^rii.i) is smaller than that for even sites, 
since doubly occupied states at odd sites may be domi- 
nated by higher order hopping processes. The variance of 
local spin fluctuations AS* shows somewhat complicated 
behavior in the intermediate U region, since it is given by 
|(( n i) — ^( n i,l n i,l))- With a large U limit, however, AS 1 ? 
is simply proportional to the local density (rij), since the 
double occupancy tends to be zero. 

We next investigate the time evolution of a band in- 
sulating state realized at half filling. Owing to the two- 
site periodic superlattice structure, the atoms are ini- 
tially distributed with a spatial configuration of 0,2,0,2, 
. . .,0 since Vd(=20) is very large. This drives the sys- 
tem to a kind of band insulating state. 17-21 To elucidate 
the dynamical properties of the band insulating state at 
half filling, we study the time evolution of the physi- 
cal quantities mentioned above. In Fig. 3, we present re- 
sults computed for several different choices of U. Let us 
first focus on the local density (rij) for even and odd 
lattice sites as shown in Fig. 3(a) and (b). The local den- 
sity at even and odd sites clearly exhibits the particle- 
hole symmetry. It is seen that the oscillation period de- 
creases as U increases. The short time region shows the 
t -dependence and its coefficient does not depend on U, 
which is also the case for quarter filling. On the other 
hand, the oscillation period in (rij) is strongly affected 
by the on-site interaction, in contrast to the quarter 
filling case. The period is approximately estimated in 
the limit of V d ^> U J, t, where the local density 
can be expressed as ^-(1 — cos(Ut)) based on time- 
dependent perturbation theory This implies that the os- 
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Fig. 3. (Color online) Similar plots of the local quantities to those 
in Fig. 2: half-filling case with L = 36 and N = 36. The on-site 
interaction is chosen as U = 1 (solid line), 5 (dashed line), and 
9 (dotted line) at (a) i = 18 and (b) i = 19. 



dilation period is changed as l/U, which is indeed in 
accordance with the numerical results. The double oc- 
cupancy (rij fjij i) shows similar properties to the lo- 
cal density as regards its oscillating properties, although 
the former exhibits ^-dependence in the short time re- 
gion. The variance of local spin fluctuations is given by 
AS? = —2{m ! -[Tii ! i)), and therefore it exhibits the 

same oscillation periodicity as for the local density and 
the double occupancy. Note that AS* exhibits similar 
profiles for even and odd sites, owing to the particle-hole 
symmetry. 

In addition to the characteristic oscillating proper- 
ties, we encounter some striking behavior at half fill- 
ing. As seen in Fig. 3, the variance of spin fluctua- 
tions AS* at i = 18 and i = 19 is suppressed as 
U increases, in contrast to the naive expectation that 
the local spin develops with increasing U. Also, all 
the quantities cease to oscillate as U increases. To 
clarify the origin of this behavior, we plot the dou- 
ble occupancy summed in the unit cell (rii,-\ni,l)sum 
(= {n i=0 dd,\n i=oddii ) + {n l=even ^n i=even ^)) in Fig. 4. It 
is notable that (ni^n^) sum (~ 0.9) is hardly changed 
with time evolution for a large U , while it shows rather 
a strong time dependence for a small U . With a sud- 
den change in Vd, two atoms occupying an odd site try 
to move apart from each other, but this process is not 
allowed for U 3> J, since the large potential energy U 
cannot be transferred to the kinetic energy. This leads 
to an unusual pairing of fermions induced by a repulsive 
interaction. The resulting bound pairs can hop around 
lattice sites via a higher tunneling process called the co- 
tunneling process. For bosons, similar bound states and 
the associated co-tunneling process were proposed and 
observed experimentally. 7 The bound pairs induced by 
the repulsive interaction are stable during time evolu- 
tion and suppress the development of local spins, thus 



reducing AS? for a large U. This naturally explains the 
unusual suppression of AS? found in Fig. 3, which is also 
confirmed by the fact that the local density profile is 
roughly given by twice the double occupancy (n^j-n^j) 
for a larger U, as seen in Fig. 3. 

Finally, to make the model more accessible to experi- 
ments, we consider the effect of the harmonic confining 
potential given by the last term in Eq.(l). Figure 5 shows 
the time evolution of the physical quantities discussed 
above. The initial ground state (t=0) is assumed to be 
in a coexisting phase including two different insulating 
regions that roughly satisfy the conditions of half- and 
quarter-filling. We have confirmed that all the quanti- 
ties in each insulating region exhibit the characteristic 
time dependence discussed above: the two-site alternat- 
ing profiles are gradually smeared after the superlatticc 
potential Vd is turned off. According to the local density 
averaged in a unit cell (ni) ave (= \{{n i=od d) + {ni= eV en)) 
shown in Fig. 5, we can see that the plateau structure is 
quite stable around the center (half- filling), while it grad- 
ually collapses around the quarter-filling regions since 
they are close to the edges of the system. Although the 
spatial extent of the half- filling region {{rii) ave = 1) is 
little changed with time evolution, its insulating charac- 
teristic gradually disappears. After the sudden potential 
change, the fermions form bound pairs and the resulting 
hard-core bosons can move via a co-tunneling process by 
keeping the number of doubly occupied sites unchanged. 
Correspondingly, the spin is not developed even in the 
long-time region, although the system is in a good con- 
dition {{riijave = 1) to realize an ordinary Mott insulator 
for t > 0. On the other hand, the quarter filling regions 
are dominated by unpaired fermions with strong repul- 
sion, which are characterized by large spin fluctuations 
and vanishing double occupancy, as seen in Fig. 5. 

It should be stressed here that the half-filling region 
specified by (rii} ave = 1, which is rather robust dur- 
ing time evolution, is spatially confined around the cen- 
ter, because hard-core bosons (paired fermions) created 
around the center experience the high pressure imposed 
by strongly correlated unpaired fermions from outside. 
This causes a unique confinement phenomenon as regards 
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Fig. 5. (Color online) Spatial distribution of the local fcrmion 
density (n^), the average fermion density {rii) ave , the variance 
of local spins AS?, and the double occupancy (n^f n^|) for (a) 
t = 0, (b) t = 1, (c) t = 2, and (d) 4 = 3. The system parameters 
are taken as U = 5.0, V c = 0.025, L = 72 and /V = 32. The 
initial state is in a coexisting phase including quarter- and half- 
filling insulating regions. Characteristic saw-toothed profiles are 
gradually smeared with time evolution. 



the paired fermions due to the strong repulsive interac- 
tion. 

In summary, by using adaptive t-DMRG simulations, 
we have uncovered several intriguing dynamical features 
of repulsive fermions in a ID optical superlattice. We 
have clarified how the correlated fermions evolve in time 
when a two-site periodic superlattice is suddenly changed 
into a normal one. For a homogeneous system at quarter- 
filling, the time dependence exhibits a similar process to 
bosons. On the other hand, in the half-filling case, more 
striking features emerge: the strong repulsive interaction 
induces fcrmion pairs, which can move around the lattice 
via a co-tunneling process. This kind of pair formation, 
which is a fermionic analog of the bosonic case, 7 in turn 
suppresses the development of local spins. We have clari- 
fied how this kind of pair hopping emerges via a crossover 
from single-particle motion in the weak- interaction limit. 
We have also addressed the effect of the confining poten- 
tial. We have shown that for an initial state consisting 
of quarter- and half-filling insulating regions, an unusual 
confinement phenomenon involving paired fermions oc- 
curs with time evolution. We expect that the nonadia- 
batic properties of correlated fermionic atoms described 
in this paper will be observed experimentally. 
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